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Abstract 

As a result of their electron structure, rare earth ions in crystals at high 
temperature emit radiation in several narrow bands rather than in a continuous 
blackbody manner. This study develops a spectral emittance model for films of rare 
earth containing materials. Although there are several possible rare earth doped high 
temperature materials, this study was confined to rare earth aluminum garnets. Good 
agreement between experimental and theoretical spectral emittances was found for 
erbium, thulium and erbium-holmium aluminum garnets. Spectral emittances of these 
films are sensitive to temperature differences across the film. Emitter efficiency is also a 
sensitive function of temperature. For thulium aluminum garnet the efficiency is 0.38 at 
1700K but only 0.19 at 1262 K. 

1 Introduction 

A selective emitter is a material that emits optical radiation in a few emission 
bands rather than in a continuous spectrum like a blackbody or a gray body (constant 
emittance). Thermophotovoltaic (TPV) energy conversion is the main application for 
selective emitters. In TPV energy conversion, the selective emitter converts thermal 
energy to near infrared radiation at wavelengths where photovoltaic energy conversion 
is efficient. A TPV system is rather simple, consisting of three main components, a heat 
source, an emitter and a photovoltaic cell array. 

The ideal selective emitter would have a single emission band with an emittance 
approaching one within the band and negligible emittance outside the emission band. 
For the photon energy or wavelength region of interest for TPV energy conversion 
(1000 to 3000 nm), an electronic transition of an atom or molecule is required to 
produce the desired radiation. However, when atoms are compressed to solid state 
densities the emission is not characterized by narrow band emission as with an isolated 
atom, but by a continuous emission spectrum. 

Fortunately, there is a group of atoms that even at solid state densities behave 
nearly like isolated atoms. These are the lanthanides or rare earth atoms. For doubly 
and triply charged ions of these elements in crystals the orbits of the valence 4f 
electrons, which account for emission and absorption, lie inside the 5s and 5p electron 
orbits. The 5s and 5p electrons “shield” the 4f valence electrons from the surrounding 
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ions in the crystal. As a result, the rare earth ions in the solid state emit in narrow bands 
rather than in a continuous gray body manner. For the temperatures of interest 
(1200 <T < 2000K) the rare earths of most interest have a strong emission band in the 
near infrared (800 < X < 3000 nm) resulting from electron transitions from the first 
excited state manifold to the ground state manifold. Because of the location of their 
emission bands, the rare earths of most interest for selective emitters are ytterbium 
(Yb), thulium (Tm), erbium (Er), holmium (Ho) and dysprosium (Dy). The spectra of the 
rare earth ions in crystals have been extensively studied. Most of this work is 
summarized in the text of Dieke. 1 

The first selective emitters investigated 2 were made by sintering rare earth oxide 
powders. These emitters showed the strong emission bands. However, emittance 
outside the emission bands was also large so that the emitter efficiency was low. In the 
late 1980’s Nelson and Parent 3 " 4 reported a large improvement in rare earth oxide 
emitters. Their emitters are constructed of bundles of small diameter (5 to 10 pm) fibers 
similar to the construction of the Welsbach mantle used in gas lanterns. The very small 
characteristic dimension of these emitters results in low emittance outside the emission 
band and thus greatly increased efficiency. The fibrous selective emitters are well suited 
to a combustion driven TPV system where the fibrous mantle surrounds the flame. 
However, for coupling to other thermal energy sources likely to be used for a space TPV 
system, such as nuclear or solar, the fibrous emitter is not so well suited. For a space 
application a planar geometry is more applicable for coupling to a nuclear or solar 
thermal energy source. As just stated, it was the small characteristic dimension that 
made the fibrous emitters efficient. Another geometry for achieving a small characteristic 
dimension and also easily coupling to any thermal source is a film. A film containing a 
rare earth on a low emittance substrate, which blocks radiation from the thermal source, 
can be easily attached to any thermal source. In addition, a film is more durable than a 
fibrous geometry. Therefore, we began a theoretical and experimental investigation of 
rare earth containing film selective emitters. 5 ' 13 

Our first attempt at producing a film selective emitter was by electron beam 
evaporation of pure rare earth oxides. However, we soon learned that film thicknesses 
on the order 1 to 10 pm were not sufficient to produce large emittance. Film thicknesses 
of 0.1 to 1 mm (100 to 1000 pm) are required. Evaporation is not applicable for films of 
that thickness so we looked for other methods. We knew that yttrium aluminum garnet 
(YAG) could be doped with rare earths and grown in single crystal form. Therefore, with 
material in this form we could cut and polish a sample to any thickness. Thus our first 
successful selective emitters were single crystal rare earth doped YAG. There are many 
other possible high temperature host materials for rare earths. Some of these possible 
hosts are shown in Table I. 

The optical properties that characterize a selective emitter are the extinction 
coefficient, a x , which is the sum of the absorption coefficient, a x , and the scattering 
coefficient, a x , and the index of refraction, n. Knowing and n, the spectral emittance, 
e x , can be calculated. In the next section the model for making this calculation will be 
developed. Figure 1 shows the extinction coefficient for single crystal erbium aluminum 
garnet (Er 3 AI 5 0 12 ).The method for obtaining will be described in the Experimental 
Measurements section. Although the results in Fig. 1 apply for erbium aluminum garnet, 
qualitatively similar results will occur for any of the possible erbium containing materials. 
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The major emission band for Er occurs at wavelength, 1500 nm, with smaller bands 
occurring at A, « 970 nm, 800 nm and 640 nm. Most all high temperature ceramic 
materials have large extinction coefficient and thus large absorptance and emittance at 
long wavelengths. For YAG and the rare earth aluminum garnets this region begins at 
X ~ 5000 nm. Obviously this long wavelength region of high emittance is undesirable for 
a selective emitter. An efficient selective emitter is one that emits most of its energy in 
the large emission band. 

In the next section the theoretical model for calculating the spectral emittance, 
e^, of a one dimensional film will be presented. Following that the experimental method 
for obtaining a x , n and e x is described. Finally, experimental and theoretical results for 
and emitter efficiency will be presented. 

2 Spectral Emittance Model 


Figure 2 is a schematic diagram of the film emittance model. Thermal energy, 
Q jn , enters through the metal substrate. Part of the thermal input leaves the film at x = d 


as radiation flux Q(d).The remaining part leaves by thermal conduction, 
As a result, the energy equation is the following. 



x=d' 


-k t h 


dT 

dx 


+ Q(x) = Q in 


( 1 ) 


Where k th is the thermal conductivity (assumed to be constant) and T is the film 
temperature. The radiation flux is defined as follows. 


Q(x) = [°°q^(x,X)d?i w/cm 2 (2) 

V O 

Where q is the radiation flux at wavelength, A, and has the units w/cm 2 nm. The 
radiation flux will always be less than the blackbody flux, o sb T s 4 , therefore defining the 
following dimensionless variables 


T 



Q = 


Q 

°sb"Ts 



results in the following energy equation. 


(3) 


dT 

dx 


+ yQ = yQin = constant 


(4) 
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Where o sb is the Stefan Boltzmann constant (5.67x1 0 -12 w/cm 2 °K) and 

y = g sb~Ts d 
k th 

Thus, if y« 1 the radiation term can be neglected and the solution to Eq. (4) is the 
following. 


T = 1 - xAT (6) 

Where, 

AT = (7) 

For most of the materials considered for rare earth selective emitters 
k th > 0.01 w/cmK, T s < 2000 K and d < 0.05 cm. Thus, y will be less than 0.2. The linear 
temperature gradient given by equation (6) will be used in the following emittance 
calculation. In a recent study 14 the complete energy equation (eq. (4)) was solved. 

Based on those results, the linear temperature variation (eq. (6)) is a good 
approximation. 

To determine the spectral emittance, e x , the radiation flux leaving the film at x = d, 
q^(d), must be calculated since is defined as follows. 

e = q^(d) = q*.(d) (8) 

®bs(^’Ts) 7ri bs (A,T s ) 


Where e bs (X,T S ) is the blackbody emissive power and T s is the substrate temperature. 


e bs(^>Ts) — ra b s(^T s ) 


27rhCo 

A 5 exp(hc 0 /AkT s )-l] 


( 9 ) 


Here i bs (A,T s ) is the blackbody intensity (w/cm 2 nm steradian), h is Plank’s constant, k is 
Boltzmann’s constant, and c 0 is the vacuum speed of light. Notice that has been 
defined in terms of the substrate temperature, T s . The spectral emittance could be 
defined in terms of the film surface temperature, T f , or some combination of T f and T s . 
However, defining e x in terms of T s means < 1 in all cases since e bs (X,T S ) > q^(d). 
This definition agrees with the usual concept of emittance. 

The following analysis is similar to that given in reference 13. To calculate q^ we 
require the radiative transfer equations for radiation intensity moving in the +x direction 
at x = d, i£ (d, cos 0), and the intensity in the -x direction at x = 0, i^ (0, cos 0). 15 
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iJ(Kd.M-) = iJ(0,M-)exp 


f K A f K d 


^ J 


+ 


S> L (K*,|i)exp 


f Kh-K* x 


M ) 


dK 




( 10 ) 


0 < pi = cos 0 < 1 


i x (0,|i) = i^(K d , 


Equations (10) and (1 
the coordinate x. 


pi)exp 


fKd] 


, P , 

J 


K r 


S> t (K*,|i)exp 


( \s * \ 


K_ 

V t* j 


dK 


M- 


( 11 ) 


1 < pi = cos 0 < 0 

1 ) have been written in terms of the optical depth, K, rather than 


K = ayx 

(12a) 

K d = a x 6 

(12b) 


Where is the extinction coefficient, assumed independent of x, and is the sum of the 
absorption coefficient, a^, and the scattering coefficient, o^. 


a X- a >i + 0 >i (13) 

Appearing in Eqs. (1 0) and (1 1 ) is the so-called source function, S X (K, pi), 15 which in the 
case of isotropic scattering is independent of p (S^(K, pi.) = S^(K)). 

The intensities i^(K d ) and i^(0) are obtained from the following boundary 
conditions at K = 0 and K = K d , assuming the boundaries are diffuse (intensities 
independent of 0). 


ix(K d ) = PfoS(Kd) at K = K d (14a) 

't(0) = pfsiii(0) + e fs i xb (X,T s ) at K = 0 (14b) 

Where p f0 is the reflectance at the vacuum-film interface, p fs is the reflectance at the 
film-substrate interface and e fs is the emittance of the substrate into the film. 

Equation (14a) states that the intensity moving in the -x direction at x = d is equal to the 
reflected intensity. At x = 0 Eq. (14b) states that i^(0) is the sum of the reflected intensity 
and radiation emitted from the substrate into the film. 

At the film-vacuum interface total reflection occurs for certain angles of incidence, 
0. At an interface between a material with an index of refraction, n j; and a material with 
index of refraction m, where n, > nj, radiation moving from i into j with an angle of 
incidence 0 > 0 M , where 0 M is given by Snell’s Law will be totally reflected. Since n f > n 0 , 
for the film-vacuum interface we have the following result for the reflectance, p f0 . 
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(15a) 


PfO = 1 f° r 0 ^ 0M 


where 


cos 2 0 M = pm = i- 



v n u 


For the case where 0 < 0 M (p > p M ) we approximate p f0 by the reflectance for normal 
incidence. 15 


Pto - 


n f ~ n 0 
n f +n 0 


0<0 M (p>p M ) 


(15b) 


Since the substrate will be a metal, n s > n f) total reflection will not occur at the 
film-substrate interface. In this case we approximate p fs by the normal reflectance for a 
metal into a dielectric. 15 


Pfs - 


(n s -n f ) 2 + n 2 , 
(n s + n f ) 2 +n 2 . 


(16) 


Where n s is the real part and n s , is the imaginary part of the substrate index of 
refraction. Since the substrate is opaque. 


e fs ~ 1 Pfs 0 7) 

Now we are prepared to calculate q^(K d ). From Eqs. (10), (11), (14a) and (14b) 
we obtain 4 simultaneous equations for the 4 fluxes q£(0), q^(0), q } +(K d ), q^(K d ). Once 
those equations are solved then q^(K d ) can be determined. 


q^K d ) = q£(Kd)-q^(K d ) ( 18 ) 

After a great deal of algebra, which is outlined below, the result for q^(K d ) and thus e x is 
obtained. 

At x = 0 the flux moving in the -x direction is the following. 


qr(0) = -2kT h (O,cos0)cos0sin0d0 = 2k\ V (0, |x)(xd(x (19a) 

A JQ=%/2 Jo 

Using Eq. (1 1 ) yields the following. 

q^0) = 2q^(K d )E 3 (K d ) + O_(K d ) (19b) 

Since the boundaries are assumed diffuse, q^(K d ) = rci^(K d ). At x = d the flux moving in 
the +x direction is the following. 
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(20a) 


qJ(K d ) = it (K d> cos 0) cos 0 si n 0d0 = 27iJgi£(K d ,^d^ 

Using Eq. (1 0) yields the following 

qt(K d ) = 2q£(0)E 3 (K d ) + O + (K d ) (20b) 

Where the diffuse boundary condition means that q4(0) = 7ii ? +(0). 

At x = 0 the flux moving in the +x direction is the following. 

/•Tt rn/2 

q£(0) = - 2^J e=jt/2 p fs i^(0, cos 0) cos 0 sin 0d0 + 2 tt J e=o e fs i bs (X,T s ) cos 0 sin 0d0 (21a) 
Using Eq. (1 1 ) for i^(0, cos 0) yields the following. 

q£(°) = Pfs[2qx(Kd)E 3 (K d ) + 0_(K d )] + e fs e bs (?i J T s ) (21b) 

The flux at x = d moving in the -x direction is the following. 
q^(K d ) = 27tJ^pf O i£(K d ,cos0)cos0sin0d0 = 2 n i£(K d ,cos0)cos0sin0d0 

7t/2 

+ Pfo f ix(K d ,cos0)cos0sin0d0 (22a) 
Using Eq. (10) for i£(K d , cos 0) yields the following 

qJ K d) = 2 q£(°) PfO E 3( K d) + ( 1 -pfo)Pm E 3 — +PfO^+( K d) + 0 _ Pfo)^M — ( 22b ) 

IPmJJ vPmv 

Equations (19b), (20b), (21b) and (22b) make up the 4 equations for the fluxes q^(0), 
q£(K d ), q£(0) and q^(K d ). Appearing in those equations are the following quantities, 

O + (K d ) = 2Kj 0 Kd S(K)E 2 (K d -K)dK (23) 

04K d ) = 27tJ o Kd S(K)E 2 (K)dK (24) 

® M f^.) = 2,t)l M J Kd S(K)E 2 (^i)dK (25) 

VPm; 0 v Pm 
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and the exponential integral E n (x) defined as follows. 


fl n-2 1 

f x) 

1 z n exp 

— 

J0 1 

V zj 


dz 


Solving the 4 simultaneous equations yields the following result. 
= q£(K d ) - q^(K d ) = 


Where 


1 2[efs e bs(^)Ts) + Pfs^-(^d)]h- + < J > +(Kd)h+ 




h_ =E 3 (K d )-ng,E 3 


OV 

V J 


\ MMy 


(26) 


h M (27) 


(28) 


h+ =1-4p,sMM E 3( K d) E 3 


V^M J 


(29) 


h M =1-4p (s E§(K d ) 


(30) 


DEN = 1-4p fs E 3 (K d ) 


PfO E 3 ( K d ) + C 1 _ PfO )m mE3 


f Lf \ 

Mm > 


(31) 


The result for q^(K d ) given by Eq. (27) is slightly different than the result in Ref. 13 
(Eq. (20)). In Ref. 13 a different substrate emittance model was used and the boundary 
condition at x = d was applied incorrectly in calculating q^(K d ). As a result, DEN given by 
Eq. (31 ) replaces D in the denominator and h M replaces D as the coefficient of <h M in 
Eq. (20) of Ref. 13. Also, in Eq. (20) of Ref. 13 e fs replaces (i\ f A\ s ) 2 (1 - p Xs ) and h + is 
given by Eq. (29) shown above. 

To proceed further with the solution for q^(K d ) and thus the source function S x 
must be determined. For no scattering the source function is given as follows. 

$x(K) = n fi b (^.T) = — e b (A,,T) (32) 

71 

Since most of the selective emitter materials we have considered are single crystal in 
structure we assume that scattering is small compared to absorption and emission. 
Therefore, for the model we assume the no scattering source function given by 
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Eq. (32) is a good approximation for the selective emitters we have considered. From 
Eq. (27) and using equation (32) in Eqs. (23) to (25) produces the following result for the 
spectral emittance when a linear temperature gradient (Eq. (6)) exists 


£a = 


q?i(Kd) _ 2nf(1-p f0 ) 
e bs (X,T s ) DEN 


^f- + 2p fs O_(K d ) 


"f 


h_+0 + (K d )h + -0 M 


v Mm j 


h M 


no scattering with linear temperature gradient 


(33) 


Where, 


0 + (K d ) 


^+(K d ) 
2rif e bs (X, T s ) 


K d ( 


e -1 


<»-(Kd) 




'1 


E 2 [K d (1-v)1 


e -1 


o ex P 

■i 


u 


1- vAT 


dv 


MM 


exp 


u 


1- vAT 


dv 


(34) 


(35) 








( u \ 

v Mm j 


2rif e bs (?i,T s ) 


= PM K d( eU 


exp 


Mm 


u 


1-vAT 


4- dv 


-1 


(36) 


Where AT is the temperature gradient given by Eq. (7) and 


u = 


hc 0 

7kT c 


(37) 


_K__ x 
K d “d 


(38) 


From Eq. (33) we see that the emittance is made up of three parts. The coefficient 
of the h_term represents the radiation leaving the substrate plus the radiation reflected 
back into the film from the film-substrate interface. This part of the emittance decreases 
with increasing optical depth, K d .The second part of the emittance, 0 + h + , represents 
the radiation emitted within the film and increases with increasing K d . The last part of e x 
is the negative term, -0 M h M , which represents the radiation with angle of incidence 
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0 > 0 M that is totally reflected back into the film at the film-vacuum interface. This part 
increases with increasing optical depth. 

The integrals in Eqs. (34) to (36) can not be carried out in closed form. However, 
in the two limiting cases of K d -> °o and K d -> 0 the integrals can be evaluated. The 
K d -» oo case approximates conditions that exist in the emission band of a rare earth 
selective emitter, while the K d -> 0 case approximates conditions that exist in the region 
between the emission band and the long wavelength cutoff region. In most cases of 
interest for rare earth selective emitters the long wavelength cutoff region begins at 
X~5 pm. Also, the substrate temperatures of interest are in the region T s < 2000 K. As a 
result, the parameter, u, (Eq. (37)) is greater than 1 . Therefore, for K d -> oo the following 
result is obtained. 13 


lim _ i 


Kh — * 


O , = - e" 


UAT 


(39a) 


lim _ i 

0_ =- 
Kd —> 00 2 


U > 1 


(39b) 


lim & - 1 ii 2 p- uAT 

if . °M - ~Z f*M e 

Kd “ ^ 00 2 


(39c) 


Also for K d -> oo, E 3 (oo) = o and therefore h + = h M = 1 and h_ = 0 so that the spectral 
emittance is the following. 


lim 
Kh - > 


= n oO _ Pfo) e 


he 

~ uAT for u = — — > 1 and no scattering 
AkT s 


(40) 


For a vacuum interface, n 0 = 1 , and negligible temperature gradient, AT = 0, the usual 
result for an opaque body is obtained (Eq. (17)). 

Now consider the case where K d = 0. In that case E 3 (0) = 1/2, 

and the spectral emittance becomes the 

^ ^ V "f J 

following. 


( 1 -Pm) = ^ 

( n ol 

\ / 2 

l n f y 




Q-pfo) 


\ 


^ Pfs 



( n > 

2 

1- 

no 

( 1_ Pfo) 


v n f J 



V n f J 


8f s U 


hc 0 

XkT c 


> 1, K d = 0 and no scattering (41) 
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In this case the emittance is determined by the substrate emittance and the reflectance 
that occurs at the vacuum-film and film-substrate interfaces. If n 0 = n f (.-.p f0 = 0) then 
Eq. (41 ) yields = e fs = e os (since n 0 = n f ) as expected. 


3. Experimental Measurements 

The optical properties required to calculate the spectral emittance, e x , with the 
emittance model are the indices of refraction and the extinction coefficient. Emissive 
performance of a selective emitter can be determined by measuring as a function of 
X. Outlined below are the experimental procedures used to determine n f , a f and e^. 


3.1 Index of Refraction and Extinction Coefficient Measurements 

Extinction coefficient and index of refraction are obtained from measured values 
of reflectance and transmittance. To cover the entire spectral region from the visible to 
approximately 1 1 ,000 nm two spectrophotometers were used, for 600 < X < 2500 nm a 
Perkin Elmer Lambda 19 and for 2500 < X < 11000 nm a Nicolet 750 FTIR. 

The one-dimensional model shown in Fig. 3 was used to relate and n f to the 
measured transmittance, T^, and reflectance, R v Assuming the reflectance, p^, at the 
two interfaces is the same then the overall reflectance, R^, transmittance, T^, and 
absorptance, A x , are the following (see Ref. 15, Ch. 19). 


Ti = 


_ q 0 4 _ taQ-Pa)' 


qp 


1 .2 2 
1- Pa x a 


(42) 



Pa 


1 + 


^aO-Pa)' 

1-P 2 AT 2 A 


Pa[ 1 + x a t a] 


(43) 


Aa=i-Ra-Ta (44) 

Where x x is the internal transmittance of the material. Assuming the material behaves 
according to Beer’s law. 


^3_ = ^ = e -^ s 

q 0 2 q 0 s 


(45) 


Where S is the path length through the material. For a beam at normal incidence angle, 
as is the case for the spectrophotometers, S = d. Using Eqs. (42) and (43) the following 
result is obtained for x x . 
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(46) 


^ = e-^ d = ±< 


(i-Rx) 


+ 4 +■ 


T? 


(i-RO 2 


Where and are the spectrophotometer measured reflectance and transmittance. 
Equation (46) was used to obtain a^from the measured values of R x , and d. The 
reflectance, p^, was obtained using Eq. (43) and x^from Eq. (46). From the index of 
refraction, n f , was obtained using the expression for normal reflectance (Eq. (15b)) and 


3.2 Emittance Measurement 

Measurement of the spectral emittance is made using the experimental 
apparatus shown in Fig. 4. An atmospheric high temperature (1700 °C) furnace is used 
to heat the sample, which is mounted on the door of the furnace. The sample has a 
platinum (Pt) foil substrate and a silicon carbide (SiC) wafer behind the Pt foil. The SiC 
wafer radiatively couples well to the furnace because of its large absorptance thus 
reducing radial temperature gradients. A thermocouple to measure the substrate 
temperature, T s , is mounted through a hole in the middle of the SiC wafer. The front 
surface temperature, T f , is measured by a thermocouple mounted on a moveable probe. 

The optical components and the front surface temperature probe are mounted on 
a moveable table. Emitted spectral intensity, i ?vS , from the sample is measured with a 
1/8 meter monochromator and three different detectors. For the visible region a Si 
detector is used. For 1 000 < X < 5000 nm an InSb detector is used and for 
5000 < X < 1 1 000 nm a HgCdTe detector is used. Also 4 different monochromator 
gratings are used to span the full wavelength range. An aperture is used to limit the area 
viewed on the sample to approximately a 4 mm diameter. The intensity leaving the 
aperture is focused on the monochromator entrance by a reflecting concentrator. 

Emitted intensity is chopped to eliminate errors resulting from stray radiation entering 
the system downstream of the chopper. Upstream of the chopper, radiation emitted from 
the hot door insulation can enter the system by being reflected off the face of the 
sample. To reduce this stray light a Pt tube, which is mounted so as not to touch the hot 
insulation, is inserted through the hole in the door and up to the sample. 

The system is calibrated by moving the optical system over to the blackbody, 
which is located exactly the same optical distance away from the monochromator as the 
sample. Thus when viewing the blackbody the current output, J^ b , of photovoltaic 
detector is the following. 


J A.b = R^Fbsibbf^Tb) ( 47 ) 

Where R x is the response (amp/watt) of the detector, F Ds is the view factor of the 
detector to the sample and i bb (X,T b ) is the blackbody intensity at temperature, T b . When 
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viewing the sample the detector output is the following. 


4s - 4^Dsixs(44) 


( 48 ) 


Combining Eqs. (47) and (48) and using the definition of spectral emittance produces 
the following result. 


_ 4(44) _ 4s *bb(^» ) 
" ibs(^Ts) " 4b ibs(^Ts) 


(49) 


Notice that e ln is different than the definition for e x given by Eq. (8). In the experiment 
the intensity \ Xs that is measured results from radiation that is emitted in the direction 
normal to the surface. Thus is really the normal spectral emittance. However, in most 
cases i^g is nearly independent of angle, which was assumed in the emittance model, so 
that q x (d) = m Xs and thus s Xn = e x . 

In deriving Eq. (49) we assumed the detector response is a function of X only 
and is independent of intensity, i^. This is a good assumption for photovoltaic detectors 
but not good for photoconductive type detectors. 


4 Emissive Property Results For Rare Earth Aluminum Garnets 

Up until the present time our efforts have been directed at rare earth doped YAG 
and rare earth aluminum garnets in single crystal form. Therefore, the experimental 
results that follow are for Er and Tm single crystal aluminum garnets. Also, results for 
garnets with both Er and Ho are presented. 

4.1 Importance of Temperature Gradient on Emissive Performance 

Since radiation is strongly dependent on temperature there will be a significant 
effect of temperature gradient, AT, on the spectral emittance. This can be shown by 
using the spectral emittance model of the previous section. Using the extinction 
coefficient data for erbium aluminum garnet, Er 3 AI 5 0 12 , shown in Fig. 1 , (Eq. (33)) 
was calculated for an emitter with a platinum substrate, a thickness d = 0.63 mm and a 
substrate temperature, T s = 1635 K. Results for the platinum index of refraction were 
obtained from Ref. 1 6. The results for 800 < X < 2000 nm are shown in Fig. 5 for 
AT = 0 and AT = 0.08. As can be seen, in the emission bands centered at 
X « 1 000 nm and X = 1 500 nm, where the extinction coefficient is large, the spectral 
emittance is greatly reduced in going from AT = 0 to AT = 0.08. Outside the emission 
bands, where the extinction coefficient is much smaller, the spectral emittance is not 
greatly effected by AT. 

Obviously AT can be reduced by decreasing the thickness, d. However, 
decreasing d also reduces the optical depth, K d = a^d, which will result in reduced e x . 
Thus varying d produces counteracting effects on e x . Decreasing d will reduce AT, which 
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will increase e x . But decreasing d will also reduce K d , which will decrease e k . As a result, 
there will be an optimum thickness, d, to obtain a maximum e^. 

4.2 Effect of Doping Level on Extinction Coefficient of Er x Y 3 _ x AI 5 0 12 

The extinction coefficient, a x , is the critical optical property for determining 
spectral emittance, e x , since the optical depth, K d = a >t d, is directly proportional to a x . 
Large will result in large e x . Therefore, to evaluate the effect of doping level on the 
emittance of erbium doped YAG we measured the extinction coefficient for a series of 
doping levels. The results of these measurements are shown in Fig. 6 for the wavelength 
range 800 < X < 2000 nm. 

As Fig. 6 shows the extinction coefficient is a monotonically increasing function of 
the Er doping level. The maximum extinction coefficients in the emission bands double 
in going from Er 1 5 Y 1 5 AI 5 0 12 to pure erbium aluminum garnet, Er 3 AI 5 0 12 . However, a 
doubling of a x in an emission band where is already large does not translate into a 
doubling of e } ,. For large and thus large K d the emittance is approximated by Eq. (40) 
which shows that is independent of a^. The benefit of a large extinction coefficient is 
that the thickness, d, can be reduced while still maintaining a large optical depth 
(K d > 1 ). Thus the smaller d will result in a smaller temperature gradient which in turn will 
increase e^, as Eq. (40) indicates. 

The effect of doping level on extinction coefficient has been established only for 
E r xY3- x AI 5 0i2. However, due to their similar atomic structure we expect the same 
increase in with increasing doping level for Ho x Y 3 . x AI 5 0 12 and Tm x Y 3 _ x AI 5 0 12 . 
Panitz 17 has found similar results forYb doping of YAG. 

4.3 Spectral Emittance of ErgAlgO^TrrigAlgO^, Er o. 3 Ho 2 . 7 AI 5 ®i 2 and 

Er 1.5 H °1.5 AI 5°12 

Figure 7 compares the measured and calculated spectral emittance of erbium 
aluminum garnet, Er 3 AI 5 0 12 , for a sample of thickness, d = 0.63 mm, with a Pt foil 
substrate. The substrate temperature was T s = 1 635 K and the measured temperature 
gradient was AT = 0.08. Figure 7(a) shows the wavelength region 600 < X < 1 0,000 nm 
while Fig. 7(b) shows an expanded view of the emission band region 800 <X< 2000 nm. 
Er has four emission bands associated with electronic transitions from the first 4 excited 
state manifolds to the ground state. The most intense band centered at X = 1500 nm 
results from transitions from the first excited state manifold ( 4 l 13/2 -> 4| i 5 / 2 )- The other 
three bands centered at 1000 nm, 800 nm and 640 nm result from transitions of the 
next th ree excited state manifo Ids ( 4 I ^ ^ 4 1 1 5 / 2 ^ 4 1 9/2 — ^ ^ 15/2 ^ E 9/2 — ^ 4 li5/2>-The 

large emittance in the region X > 5000 nm results from the vibrational modes of the 
crystal lattice. As pointed out in the Introduction, this long wavelength cutoff region 
occurs for most all high temperature ceramic materials. 

There is good agreement between the measured and calculated (Eq. (33)) 
spectral emittance except for the region 1 800 < X < 6000 nm. In the region 
1800 < X < 4000 nm where the extinction coefficient is small the main contribution to the 
emittance comes from the Pt substrate. Error in the measured extinction coefficient is 
the reason the calculated is larger than the measured value in the 
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1800 < X < 4000 nm region. In this region where Er 3 AI 5 0 12 absorption, A x =1 -T^- R }l , 
is small means that T } * 1 - R v Thus, referring to Eq. (46) we see that if * 1 - R x the 
quantity T^ 2 - (1- R^r can have considerable numerical error which then results in a x 
being overestimated in the 1 800 < X < 4000 nm region. 

For the region 4000 < X < 6000 nm the vibrational modes of the crystal lattice are 
the primary source of the emission. The extinction coefficient was measured at room 
temperature. Therefore, the calculated for 4000 < X < 6000 nm corresponds to room 
temperature conditions. Since the calculated is less than the experimental in the 
4000 < X < 6000 nm region we conclude that the crystal lattice structure is changing 
with temperature such that higher energy (shorter wavelength) modes exist. These new 
modes produce the increase in extinction coefficient in the 4000 < X < 6000 nm region. 

In the emission band region, 600 < X < 1 800 nm electronic transitions of the Er ions 
account for the radiation. Therefore, widening of the emission bands is the expected 
effect for increasing temperature rather than increasing extinction coefficient. This 
conclusion is substantiated since there is good agreement between experimental and 
calculated for 600 < X < 1 800 nm with the experimental emission bands being 
somewhat broader. The resolution of the extinction coefficient data was much greater 
than experimental data. As a result, the calculated e x shows more structure in the 
emission bands. 

As a measure of the effectiveness of Er 3 AI 5 0 12 as a selective emitter we define 
the emitter efficiency as follows. 


useful radiated power _ j 0 c lx( d ) cl ^ _ J Q e^bsiXTs) 1 ^ 
total radiated power = J o " qx(d)dx . = J o “ e ,e bs (X,T s )ca 


The numerator is the power radiated in the wavelength region 0 < X < X € , where X € is the 
wavelength at the end of the main emission band (X € = 1 650 nm for Er 3 AI 5 0 12 ). The 
denominator is the total radiated power. As shown in Fig. 7 the calculated efficiency is 
t|e l ca | = 0.20 and the experimental efficiency is ri E | exp = 0.21 . The theoretical efficiency 
is smaller because e^| ca | is too large in the 1 800 < X < 4000 nm region. As Eq. (33) 
shows, is not a function of T s but is a function of AT for the no scattering and linear 
temperature gradient approximations. Flowever, the blackbody emissive power, 
e bs (X,T S ), is a sensitive function of T s . As a result r| E is also a sensitive function of T s . 
For d = 0.3 mm, T s = 1234 °K and AT = 0.094 the measured efficiency is reduced to 

T lElexp = 0 - 065 - 

Figure 8 compares calculated and measured values of e^forTm 3 AI 5 0 12 on a Pt 
substrate with d = 0.5 mm. The substrate temperature wasT s = 1700 K and AT = 0.1 1 . 
Similar to Er 3 AI 5 0 12 , Tm 3 AI 5 0 12 has emission bands that originate from transitions to 
the ground state from the 4 lowest excited states; 3 F 4 3 PI 6 centered at X- 1700 nm, 
3 H 5 ^ 3 PI 6 centered at x= 1200 nm, 3 PI 4 -^ 3 H 6 centered at X= 800 nm and 3 H 3 -^ 3 H 6 
centered at X = 700 nm. The fi rst excited state to ground state transition ( 3 F 4 — > 3 H 6 ) has 
the largest intensity. Also similar to Er 3 AI 5 0 12 , the long wavelength region of large 
emittance begins at X = 5000 nm. Similar to Er 3 AI 5 0 12 , there is good agreement 
between the experimental and theoretical emittance in the emission band regions. Also, 
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because of experimental error in the extinction coefficient for the region 
2000 <X< 5000 nm the theoretical emittance has large fluctuations, as well as, being 
too large. ForTm 3 AI 5 0 12 the measured emitter efficiency using X e = 1900 nm in Eq. (50) 
is r| E | exp = 0.38 compared to r| E | ca! = 0.31 for the sample in Fig. 8. Since T s is nearly the 
same for both the Er 3 AI 5 0 12 and Tm 3 AI 5 0 12 sample we can compare their efficiencies. 
The larger rj E forTm 3 AI 5 0 12 compared to r| E for Er 3 AI 5 0 12 results for two reasons. First 
of all, even though in the four emission bands is nearly the same for both Tm and Er 
the widths of the Tm bands are greater. Second, the integration limit, X ( , in Eq. (50) is 
larger forTm (X € = 1900 nm) than Er = 1650 nm). AtT s = 1262 K and AT = 0.062 the 
measured efficiency is r| E | exp = 0.19 for the sample in Fig. 8. 

One possible method for increasing the useful power (numerator of Eq. (50)) is by 
doping the host material with more than one rare earth. To investigate this approach we 
measured the spectral emittance and efficiency of erbium holmium aluminum garnet 
(Er x Ho 3 . x AI 5 0 12 ) at two different doping levels; Er 03 Flo 27 AI 5 0 12 and Er^FIOj 5 AI 5 0 12 . 

The experimental and theoretical results for Er 03 hio 27 AI 5 O 12 of thickness, d = 0.55 mm, 
with Pt foil substrate are shown in Fig. 9. In this case the substrate temperature was 
T s = 1 298 K and the temperature gradient was AT = 0.080. Notice there is good 
agreement between measured and calculated e^. Since the doping level of Er was much 
lower than Flo the four emission bands of Er at i = 1 500, 1 000, 800, and 640 nm are 
smaller than the Flo emission bands resulting from transitions from the four lowest 
excited states. The main Flo emission band results from transitions from the first excited 
state manifold to ground state, 5 I 7 5 I 8 , and is centered at 2000 nm. Transitions 

from the second excited state manifold, 5 I 6 5 I 8 , produce the band centered at 

X * 1 1 00 nm. Transitions from the third excited state manifold, 5 I 5 -> 5 I 8 , produce the 
band centered at 890 nm. And finally, transitions from the fourth excited state 
manifold, 5 I 4 5 I 8 , produce the band at X = 750 nm. The three shorter wavelength bands 

of Er and Flo are close together so that they merge together in the spectral emittance. 
Using X t = 2150 nm in Eq. (50), the experimental and calculated efficiencies are 0.27 for 
the data in Fig. 9. 

When the doping levels of Er and Flo are the same (Er, 5 Flo 1 5 AI 5 0 12 ) the 
spectral emittance results are those shown in Fig. 1 0 for a sample of thickness, 
d = 0.65 mm, and Pt foil substrate. In this case, the substrate temperature was 
T s = 1 21 8 K and the temperature gradient was AT = 0.074. The Er bands have 
increased e x while the Flo bands have decreased e x compared to the Er 0 3 FHo 2 7 AI 5 0 12 
results (Fig. 9). Again there is good agreement between e^| ca[ and ejJ exp . The measured 
efficiency is r| E | exp = 0.23 compared to rj E | ca | = 0.20. Since T s and AT are nearly the 
same for Er 03 Flo 2 7 AI 5 0 12 and Er., 5 Flo 1 5 AI 5 0 12 in Figs. 9 and 10, we can make a valid 
comparison of efficiencies. Since r| E is larger for Er 0 3 Flo 2 7 AI 5 0 12 , we conclude that it is 
a better selective emitter than Er 1 5 FHo 1 5 AI 5 0 12 . 


5 Conclusion 

Rare earth ions in high temperature host materials emit in several narrow bands 
in the near infrared and the visible spectrum. A spectral emittance model, where optical 
depth is the key variable, has been developed for these rare earth selective emitters. 
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There are many possible host materials although we have investigated only rare earth 
aluminum garnets of Er, Ho and Tm. Experimental measurements of the extinction 
coefficients and spectral emittances of these garnets have been made. The 
experimental and theoretical spectral emittance results are in good agreement. 
Maximum spectral emittances of 0.6 to 0.8 occur in the emission bands. At X > 5000 nm 
all these materials have large emittance (e 3l > 0.8) as a result of emission from the 
vibrational modes of the crystal structures. Small temperature gradients (AT < 0.1) 
across the film emitters result in significant reductions in spectral emittance. Because 
spectral emittance increases with optical depth but decreases with temperature gradient 
there is an optimum film thickness for maximum emittance. Emitter efficiency is a 
sensitive function of temperature. For thulium aluminum garnet (Tm 3 AI 5 0 12 ) the 
efficiency atT s = 1262 K for film thickness, d = 0.5 mm, and temperature gradient, 

AT = 0.062 is r| E | exp = 0.19 and atT s = 1700 K with AT = 0.1 1 the efficiency is 
Ti E lexp= 0-38- 
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TABLE I.— POSSIBLE HIGH TEMPERATURE HOST 
MATERIALS FOR RARE EARTHS (Re) 


Host material 

Chemical formula 

Melting point, °C 
of host 
(Ref. 16) 

Pure rare 
earth oxide 

Re 2 0 3 

>2200 

Yttrium aluminum 
garnet (YAG) 

R e xY3-x^50i2 

0<x<3 

1930 

Gadolinium gallium 
garnet (GGG) 

Re x Gd 3 .xGa 5 0i 2 

0<x<3 

1750 

Spinel 

Re x Mgi_xAI 2 0 4 

0<x<1 

-2100 

Zirconia 

RexZr^O, 

0<x<1 

-2700 

Yttria 

^®x^2-x^3 

0<x<2 

2400 
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Fig. 1 . Extinction coefficient for erbium aluminum 
garnet (E^AbjO^)- 


Vacuum, n Ao = 1 



p fo = reflectance at film - vacuum interface 
p fs = reflectance at film - substrate interface 
£ fs = emittance of substrate into film = I - p fs 
i bs ( A, T s ) = blackbody intensity for T = T s 

Fig. 2. Schematic diagram of film emittance model. 
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Fig. 3. Model for calculating overall reflectance, 
R\, and transmittance, T^. 
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Fig. 4. Experimental emittance measurement 
apparatus. 



Wavelength, X nm 

Fig. 5. Effect of temperature gradient on spectral 
emittance of Er^At^O^ with platinum substrate for 
T s = 1635 K and d = 0.63 mm. 
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Fig. 6. Effect on extinction coefficient of El -doping 
level in YAG. (a) Extinction coefficient for 
Ei'i 5Y 1 5AI5O [2- (b) Extinction coefficient for 
Ei'i sY 12 A I5 0 1 2 - (c) Extinction coefficient for 
EnjYojAI^O 1 2. (d) Extinction coefficient for 

Er^AI^O^. 
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Wavelength, X, nm 



Wavelength, X, nm 


Fig. 7. Comparison of experimental and calculated 
spectral emittance of Erbium aluminum garnet 
(Er^Al^Ojj) with platinum substrate. Thickness, 
d = 0.63 mm, substrate temperature, T s = 1635 K, 
front face temperature, Tf =1503 K, AT = 0.0807. 
(a) 600 < A. < 10,000 nm. (b) 800 < A < 2000 nm. 
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Wavelength, X, nm 


Fig. 8. Comparison of experimental and calculated 
spectral emittance of Thulium aluminum garnet 
(T 1 TI 3 AI 5 O 12 ) with platinum substrate. Thickness, 
d = 0.50 mm, substrate temperature, T s = 1700 K, 
front face temperature, Tf =1517 K, AT = 0.0108. 
(a) 600 < A. < 10,000 nm. (b) 650 < \ < 2500 nm. 
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Wavelength, X, nm 



Wavelength, X, nm 


Fig. 9. Comparison of experimental and calculated 
spectral emittance of Ero. 3 Ho 2 , 7 Al 50 i 2 with 
platinum substrate. Thickness, d = 0.55 mm, 
substrate temperature, T s = 1298 K, front face 
temperature, Tf = 1 194 K, AT = 0.08. (a) 600 
< A. < 10,000 nm. (b) 800 <\< 2500 nm. 
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Fig. 10. Comparison of experimental and calculated 
spectral emittance of Er ] 5 H 01 5AI5O12 with 
platinum substrate. Thickness, d = 0.65 mm, 
substrate temperature, T s = 1218 K, front face 
temperature, Tf =1 128 K, AT = 0.074. (a) 600 
< A. < 10,000 nm. (b) 600 < A < 2500 nm. 
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